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Abstract 

We study the localization length l c of a pair of two attractively bound particles 
moving in a one-dimensional random potential. We show in which way it depends 
on the interaction potential between the constituents of this composite particle. 
For a pair with many bound states iV the localization length is proportional to N, 
independently of the form of the two particle interaction. For the case of two bound 
states we present an exact solution for the corresponding Fokker-Planck equation 
and demonstrate that l c depends sensitively on the shape of the interaction potential 
and the symmetry of the bound state wave functions. 
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1 Introduction 



Within the single-parameter scaling hypothesis [1] the dimensionless conduc- 
tance g(L) as a function of system size L determines the transport properties of 
noninteracting electrons in a random potential (for a review, see Kramer and 
MacKinnon [2]). In one-dimensional disordered systems all states are local- 
ized [3] and the conductance g(L) decreases exponentially [4] with increasing 
sample size L. The localization length l c can be derived from the asymptotic 
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behavior of the conductance g(L) using the relation 

l? = - lim ^-lng(L). (1) 

The localization length l c and the Ljapunov exponent 7 ~ Z^ 1 are self-averaging, 
i.e. non-random quantities [5] whereas the distribution function of the conduc- 
tance is approximately log normal [6]. In the weak disorder limit characterized 
by a small Fermi wavelength compared to the scattering length I the local- 
ization length l c is given by the mean free path / for backward scattering 
[5]. 

In this paper we study the localization of a composite particle with internal 
degrees of freedom. As a model we consider a single pair of two particles in 
a weak one-dimensional white noise potential. We restrict our considerations 
to the limit of two tightly bound constituents so that the pair does not decay 
during its motion through the disorder potential. Considering the semiclassical 
limit the kinetic energy of the center of mass motion is assumed to be larger 
than the typical energy of the disorder potential. Typical length scales of our 
model include the mean free paths lb and // of the pair for backward and 
forward scattering, respectively. The scattering processes cause a transfer of 
kinetic energy of the center of mass motion to the relative motion within the 
pair and vice versa. Therefore, the scattering lengths If, and If depend on the 
structure of the pair, i.e. on the interaction potential between the constituents. 
Thus one expects that the pair localization length l c also depends on the shape 
of the attractive pair potential. 

As another but different example for the interplay between interaction and 
disorder many recent papers study the coherent propagation of two interact- 
ing particles (TIP) in a one-dimensional random potential (see [9] and the 
literature cited there). It is worth mentioning that the TIP model considers 
the propagation of two repulsing or weakly attracting particles which form a 
pair due to the scattering by the disordered potential. The size of the pair 
is of the order of the one particle localization length /. As Shepelyansky [10] 
and Imry [11] argued, attractive as well as repulsive interaction between the 
particles cause a localization length of the pair which is much larger than the 
single particle localization length. In the present work, we study two parti- 
cles with strong interaction which form bound states and have a size that is 
much smaller than the one particle localization length. For that reason the 
TIP model and the composite particle model studied in this work consider 
different limiting cases, which are not directly comparable. Moreover, in the 
mentioned TIP approach a lattice model of two particles interacting at one 
lattice site is considered. Generally the coherent motion of the interacting pair 
is studied only in the middle of the band (E — 0). However, we study a con- 
tinuum model of two interacting particles in the semiclassical regime with a 
Fermi wavelength much smaller than the scattering length (\ F <C I). 
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A composite particle with iV bound states may be considered as an iV-channel 
problem for the center of mass motion. On the other hand the motion of a sin- 
gle electron in a thick wire is also a multichannel problem since the transversal 
motion is quantized. In the thick wire the scattering processes between dif- 
ferent channels are random and of the same strength. The localization length 
l c of the wire increases with the channel number A as l c ~ Nl because effec- 
tively only one of these channels gives rise to coherent backward scattering 
and produces localization [7,8]. In our composite particle model the scattering 
processes between different bound states are also random. However, in con- 
trast to the model of a thick wire, the scattering probabilities depend strongly 
on the channel numbers. Therefore it is interesting to investigate the localiza- 
tion length in the composite particle model and compare the results for l c of 
both different models. 

The model of two particles with strong attractive interaction might be relevant 
for the Coulomb-correlated electron-hole pairs in disordered semiconductors. 
Recent work [12] has dealt with the dynamics of the electron-hole pair in a 
ID model. Another physical problem which can be related to the model we 
consider below is the superconductor-isolator transition. In order to study the 
effects of interaction and disorder, Lages and Shepelyansky [13] investigated 
numerically the Cooper problem of two quasiparticles with attractive interac- 
tion above a frozen Fermi sea in the presence of disorder. 

The basic ideas of the model we consider were introduced by Dorokhov [14]. 
He studied a pair of particles bound by a harmonic oscillator potential in a 
weak random disorder potential. For this very special interaction he found 
a dominant forward scattering in the case of many relevant bound states A 
{h/ If ~ In A) employing the specific properties of the harmonic oscillator 
eigenfunctions. Moreover, the pair scattering length l^ was shown to be on the 
order of the single particle mean free path I and independent of the strength 
of the oscillator potential. Using the fact that the forward scattering is the 
dominant process the pair localization length l c = Nlb/2 ~ Nl was calculated. 

Our results show that Dorokhov's conclusions are generic for any attractive 
interaction potential in the case of many bound states A 3> 1 regardless of its 
shape. We extend the method [14] from the harmonic oscillator potential to 
an arbitrary attractive interaction potential. This can be accomplished since 
the excited states of the relative motion within the pair can be described in 
the semiclassical approximation. We show that the specific properties of the 
harmonic oscillator wave functions (in [14] the result was obtained using the 
recurrence relations of the Hermitian polynomials) are of no importance in the 
multi-channel case. Using WKB wave functions for the bound states of the pair 
we estimate the scattering matrix elements and the mean free paths lb and If 
in the case A ^> 1. We find that lb^$> If for any interaction potential. We also 
show that 1^ is qualitatively independent of the interaction. For the localization 
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length l c we obtain l c ~ Nib/ 2 which implies that l c is also independent of the 
interaction. 

In addition to the considerations for many bound states iV > 1 we also 
investigate the first nontrivial case of a small number of bound states, namely 
N = 2. For a single particle in a random white noise potential the mean 
free paths lb and If are equal. We show that already for a pair with only 
two bound states, the forward scattering is enhanced relative to the backward 
scattering and that the ratio If /lb depends sensitively on the energy. Finally 
we derive an exact expression for the localization length of the pair with two 
accessible bound states. In contrast to earlier treatments [7] of this problem 
we do not have to impose any additional restrictions to the scattering matrix 
elements. Therefore, our solution allows to study the influcence of the two 
particle interaction on the localization directly. The result indicates that no 
simple relation exists between l c and the mean free path of the pair l b as it 
was in the limit of many channels N ^> 1. As an application of our general 
result we present numerical calculations showing that the localization length 
l c depends sensitively on the interaction potential, the symmetry of the bound 
state wave functions and the total energy of the pair. 

The paper is organized as follows. In Section 2 we describe the model and 
briefly summarize the transfer matrix approach applied to the composite par- 
ticle model with arbitrary interaction potential between the two constituents. 
In Section 3 we consider a pair with high enough total energy allowing for 
many bound states. We determine the scattering lengths l b and If as well as 
the pair localization length l c . In Section 4 we derive an exact expression for 
the localization length of a pair with two bound states (N = 2). Furthermore, 
we present numerical results for the localization length l c as well as for the 
mean free paths lb and If for different interaction potentials. 



2 Model and method of calculations 

We consider two interacting particles in a one-dimensional random potential 
V(x). The Hamiltonian is given by 

H = -\^-^-2+ <V) + V(* + y/2) + V(x - y/2), (2) 

where the units are chosen in a way that H — 1 and the single electron mass 
m = l. The variable x denotes the center of mass coordinate of the pair while 
y represents the relative coordinate. The symmetric interaction potential u(y) 
between the two particles is assumed to be attractive. The bound states (f> n (y) 
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and the corresponding energy eigenvalues e n are given by 

d 2 



+ u(y) 



^>n(y) = e n (f) n (y). 



(3) 



In the case of two identical particles the states 4> n {y) of even parity describe 
two bosons or a singlet state of fermions while the wave functions of odd 
parity correspond to a triplet state of fermions. The random potential V(x) 
is nonzero in the finite interval < x < L. This region of length L is enclosed 
by two ideal leads without any disorder. 

The scattering of the center of mass of the pair with N bound states can be 
mapped onto an iV-channel problem of a single free particle. The stationary 
scattering states of the pair in the leads are given by 



7V-1 

*(x,y) = ]T My) 



n=0 



A 



exp (ik n x) 

\/kr, 



+ B, 



exp (—ik n x) 
\/k n 



(4) 



Here, the momentum of the center of mass is denoted by k n . The total energy 
of the pair E = k^/A + e n determines the number of bound states (channels) 
iV which are included in the sum in Eq. (4). For a given total energy E this 
number iV can be obtained from the relation e^r-i < E < e^. 

As far as the coefficients A n and B n are concerned one has to distinguish be- 
tween the left lead and the right lead. They can be written as vectors (A L , B L ) 
and (A R , B R ) for the left and right one, respectively. The transfer matrix r 
describes the scattering in the region between the two leads. It relates the 
states in the right lead to the states in the left lead [8]: 



(5) 



The transfer matrix r is a 2N x 2N symplectic matrix which can be expressed 
by the NxN reflection matrix r and the transmission matrix t. These matrices 
describe the scattering between the channels n and m which correspond to the 
associated bound states of the pair. The dimensionless conductance g can be 
expressed [15] in terms of the N eigenvalues T n of the Hermitian matrix t't as 
g — T n . Therefore, it is convenient to employ the Hermitian matrix M = rV 
which can be expressed by the polar decomposition as [7] 





M = 



u 
u* 



cosh2r -sinh2r 
-sinh2r cosh2r 







It 


u 1 ) 



(6) 



In this parameterization u is a iV x iV unitary matrix and r is a iV x iV diag- 
onal matrix. The transmission eigenvalues T n are then related to the matrix 
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elements T n by T n = 2/[cosh(2r n ) + 1]. 

To evaluate the transmission eigenvalues T n (L) for a given size L we now 
consider the matrix M as a function of the sample size L. As was shown 
in [14] the evolution of M(L) with increasing L can be described by a set 
of Langevin equations. The random potential V(L) enters these equations 
as a multiplicative factor. In order to derive these equations the following 
assumptions were made: 

• The potential V(x) is a weak Gaussian white noise potential with the cor- 
relation function 

(V(x)V(x')} = 2D5(x - x'), (7) 

where (...) denotes the ensemble average. The mean free path I of the pair 
in comparison to the inverse wave numbers k~ l is the appropriate quantity 
to describe the strength of the disorder. The semiclassical limit is considered 
where k n l ^> 1. The mean free path / of a single particle with Fermi energy 
E F is given by I — E F /D. 

• The mean free path I is large compared to the typical width b n of the bound 
states <j) n (y): I b n . This assumption ensures the concept of mean free path 
for the center of mass motion. 

The localization length l c of the pair is studied by the asymptotic behavior of 
the transmission eigenvalues T n (L) for large system sizes L — > oo. These coeffi- 
cients T n (L) decrease exponentially with the system size L. The corresponding 
length scales are defined by 

l? = - lim -^(lnT n >. (8) 

The transmission eigenvalues can be arranged in a hierarchical order T/v-i 
Tn-2 C ••• CT € 1. The localization length l c is then given by the largest 
length l . All lengths l n averaged over the channels just give the mean free 
path l b of the center of mass for backward scattering [7]: 

-i N-l 

j; E C 1 = K 1 0) 

iV n=0 

(see Eq. (14) below for the backward scattering length). Taking into account 
the exponential smallness of all transmission eigenvalues and their different 
order of magnitude, the Langevin equations for the matrix elements T n and 
u nm are given by 
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Br 

-Q^ = ii^-^) m -V(L) (10) 
du N ~ 1 

W Ll ri m 



2i J2 [ a mui m + u nl (P*)i m e(l - m) 



+ u nl ((3 u ) lm Q(m-l)]-V(L) (11) 

where the abbreviation (3 U stands for (3 U = u^(3u*. The value of the step func- 
tion Q(k) for k = is defined as 0(0) = 1/2. The matrix elements a nm and 
f3 nm in Eq. (11) describe the forward scattering and the backward scattering 
of the pair, respectively. They depend on the two particle interaction and can 
be written as 



&nm 1 11 ~ ( ^ n km) ■> 



h m 



Pnm = j , F W nm {k n + k m ), (12) 

where W nm {k) is 

W nm (A;) = 2 J dy <t>n{y)Mv) cos . (13) 

Within the Born approximation, the mean free paths l^j of the pair for back- 
ward and forward scattering are given by [14]: 

gjj N-l gjj N-l 

l b l = jE \Pnm\ 2 and Ij 1 = — ]T Km\ 2 - (14) 

n,m n,m 

Note that one has to consider bound states of a given parity only. Due to 
the finite width of the pair we therefore expect that the forward scattering 
is enhanced relative to the backward scattering although the scattering of a 
single particle is isotropic. This circumstance allows to solve the many channel 
problem [14] with iV 3> 1 (see section 3). 

Using the fact that the transport is effectively determined by just one channel 
(T ^> T n>0 ) Eq. (8) results in the following expression for the pair localization 
length [14] 

n-i y 
Z" 1 = 16D (Ko| 2 Ko| 2 >(l - -S nm )\Pnm\ 2 - (15) 

n,m=0 

Hence, the localization length l c depends on the correlation function (|w„o| 2 |w m o| 2 ) 
which has to be evaluated using the nonlinear Langevin Eqs. (11). 



7 



For N — 1 the unitary matrix u is simply a phase factor with \uqq\ — 1- 
The localization length / c is then given by the mean free path l} } for backward 
scattering 

C 1 = K 1 = 8L>|/3oo| 2 (16) 

and the interaction potential «(?/) enters via the ground state wave function 
in the matrixelement (13). In the next section we consider the multi-channel 
case iV 3> 1 and show that the localization length depends only weakly on the 
interaction. 



3 Multi-channel localization length 



In order to evaluate the localization length (15) in the multi-channel case 
N ^> 1 one has to compute the correlation function (|w n o| 2 |wmo| 2 ) using the 
nonlinear Langevin equations (11). For the harmonic oscillator interaction it 
was shown in [14] that the forward scattering dominates in Eq. (11). Therefore 
the nonlinear terms due to the backward scattering processes were neglected 
in (11) and it was found that the correlation function is properly described by 
the invariant unitary ensemble 

(KolVrnol 2 ) = (l+5 nm )/N 2 . (17) 

In the following we proof that the same result is valid for any attractive inter- 
action potential. The order of magnitude of forward and backward scattering 
processes is determined by the corresponding mean free paths If and given 
by (14). These mean free paths depend on the pair interaction potential u(y) 
via the bound state wave functions 4> n {y) in Eq. (13). We calculate the ma- 
trix elements W nm (k) using WKB wave functions and evaluate the integral 
(13) in the saddle point approximation. This approximation is justified by the 
semiclassical approximation and the assumption N 1, which implies that 
the vast majority of bound states is properly described by the WKB approx- 
imation. In this approximation the following expression can be derived from 
(13) 



W nm {k) = 4 (-1)(»±"0 



2tt 



1/2 



T n T m \u'(y ) ■ k\ 



x cos 



\k\ 



7Y 



T- / dky {k)-- 



&2,1 



(18) 



The time T n denotes the period of the semiclassical motion of the bound state 
n. The two saddle points ±y are functions of k and follow from the implicite 
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Fig. 1. Comparison of the exact function W22(k) and the calculated WKB estimate 
in the case of the harmonic oscillator interaction. The WKB approximation is ex- 
cellent already for small quantum numbers in the wave function (in this example 
4>n= 2 ) as l° n g as k lies in the oscillating regime. The WKB approximation clearly 
breaks down at = and k 2 ^ = 12. 

equation 

/ \ k (e ra Cm) /in\ 

«fa>) = " is " —0— ( 19 ) 

In Eq. (18) u'(y ) stands for the derivative of the potential evaluated at po- 
sition i/q. It results from the Gaussian integral at the saddle point. Calcu- 
lating the matrix elements W nm (k) for backward or forward scattering we 
set k = k n ± k m . In this relation as well as in Eq. (18) the upper sign and 
k 2 characterize the backward scattering while the lower sign and k\ repre- 
sent the forward scattering. The expression (18) is valid within the interval 
k\ < \k\ < k 2 where Eq. (19) has a real solution y (k). The limits of the 
interval are given by 

k 2 ,i = 2(\/en - u ± Vem - uo) (20) 

with the potential minimum uq = u(y = 0). Outside this k region the saddle 
point is shifted to the complex y-plane and the matrix element W nm (k) de- 
creases exponentially [17]. Calculating the mean free paths (14) we will neglect 
these exponentially small contributions. 

To demonstrate the accuracy of the WKB wave functions together with the 
saddle point approximation we show an example for a matrix element in Fig. 
1. We have calculated the exact function W 2 2(k) for the harmonic oscillator 
potential u osc (y) = £,~ 2 (y/£,) 2 using the analytical expression for the third 
symmetric bound state 4> 2 (y). This exact result is compared with our WKB 
result (18). One clearly sees that the oscillatory part is very well described by 
the approximation. Of course, close to k — ki t2 the saddle point approximation 
fails because u'(y ) goes to zero and the Gaussian approximation is not valid. 
Furthermore, the diagonal matrix elements W nn cannot be calculated from Eq. 
(18) for k — 0. However, definition (13) immediately yields W nn {k = 0) = 2. 

The scattering lengths are then evaluated by replacing the sum in (14) by an 
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integral (with k n — > k and k m — > k). We obtain 

r 1 ^ [ ^ ^ (21) 
bJ = 4^NJ (k±k)- \u'{y )\' ( ' 

This result may be understood as follows. For k 2 = (k n ±k m ) 2 the saddle point 
equation (19) can be rewritten as 

k 2 P P P 

u{ m )=E-^- k -f = e n - k -f = e m -^. (22) 

This means that the kinetic energy of the relative motion e n — u(yo) for the 
initial channel n transforms into the kinetic energy P^/4 of the center of 
mass in the final channel m and vice versa. The transition rate into the final 
channel m is therefore proportional to the probability density p n {k m ) of finding 
the pair in the initial channel n with the appropriate relative momentum. 
The expression (21) for the mean free paths follows then immediately from 

Pn(k m ) =4/(T>'(y )|)- 

We will now show that Eq. (21) implies a mean free path for the backward 
scattering that is almost independent of the channel number N. Using Eq. 
(22) we write u'(y ) formally as 



u'(y ) = F(u(y )) — F [E - — . (23) 




This is always a unique representation if the interaction potential is assumed 
to be monotonic. Introducing polar coordinates the angular integration in (21) 
can be performed and we obtain the expression for the backward scattering 
length 

i 

~V2 + l 



IZ 1 = V2\n 



!dt{i-t)-^ [FiEtT 1 

o 



-n 1 . (24) 



jdt(i-ty/* [F(Et)]- 1 



The denominator in (24) results from the number N of channels with given 
parity below the energy E. Within the WKB approximation the relation be- 
tween the number iV and the energy E can be expressed by using expression 
(23) once more: 

N(E) = ±-jdy^E- u(y) 
1 )■ 

= -E^ 2 dt(l- t 2 ) 1 ' 2 [F(Et)]- 1 . (25) 

7T J 



From Eq. (24) follows that If, ~ I with a prefactor of order one. This prefactor 
depends only weakly on the interaction potential u(y) and the number of 
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accessible bound states N. To demonstrate this we consider two examples of 
interaction potentials. 



In the case of a scale invariant potential with u(y) ~ y v and F(u) ~ u^ u ^ l ^ u 
one immediately realizes that lb /I is completely independent of the energy as 
it was the case for the harmonic oscillator [y = 2). From (24) follows 



r 1 = rHV2\n 



V2 + 1 
y/2- 1 



5 + -). (26) 



This expression includes the result for the harmonic oscillator presented in 
[14] which was found by calculating the matrix elements W nm exactly. For 
any interaction potential u(y) with many bound states we expect only a weak 
dependence of 4 on the energy E or on the number of channels N which are 
below this pair breaking energy. As a further example we give the result for 
the Poschel- Teller interaction [16] u a PT (y) = £ _2 (a 2 — a(a + 1)/ cosh 2 [y/£]) 
where a is an integer parameter. In contrast to the scale invariant potentials 
this interaction has a maximal number N max of bound states given by the 
parameter a. The pair breaking energy is ~ £~ 2 a 2 . The semiclassical limit 
can be achieved for a large parameter a implying N max ^> 1. The backward 
scattering length (24) can then be written as 



1- 1 = r^in 



V2 + 1 
V2-1 



(N/N max ) 



2 ' 



(27) 



If the total energy of the pair is much smaller than the pair breaking energy 
then iV N max follows and again the ratio lb/ 1 is only very weakly energy 
dependent. 

Using the same approach for the forward scattering length If one realizes that 
the main contribution in Eq. (21) is due to the scattering between adjacent 
channels with k ~ k. However, the saddle point approximation fails for k = 
k and the integral in (21) diverges. Therefore, one has to exclude a small 
strip of order \k — k\ ~ 1/N in the integral and calculate the in-channel 
scattering separately. We find in accordance with the ealier result for the 
harmonic oscillator interaction [14] that 

Ij 1 ~ r 1 In iV (28) 

for iV ^> 1. From this estimate (28) follows If <C lb for the multi-channel case. 
This means that the forward scattering is dominant and the relation (17) can 
be applied for any interaction potential u(y). The Eqs. (24, 28) are thus the 
first major result of this work. 
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Using the correlation function (17) we then obtain the localization length 

N 

lc = -jh (29) 

from expression (15). Thus we find a delocalization proportional to the num- 
ber of open channels N compared to the one particle localization length which 
equals I. This shows that the localization effect is qualitatively independent 
of the shape of the attractive interaction potential u(y). The interaction only 
gives rise to a weak energy dependence and a numerical prefactor in the back- 
ward scattering length of the pair l b compared to the single particle scattering 
length /. 

Our result (29) has the same structure as in the different model of a thick 
wire [8] . For that case the Dorokhov-Mello-Pereyra-Kumar (DMPK) equation 
describes the evolution of the distribution function for the transmission eigen- 
values T n (or r„) with increasing length of the wire. Compared to our model 
an important simplification there is the isotropy assumption for the scattering 
between the channels on average. Using this assumption the Fokker-Planck 
equation for the transmission eigenvalues T n decouples from the distribution 
function of the matrix elements u nm . At least in the absence of time-reversal 
symmetry the DMPK equation could be solved exactly [18] . In our model the 
dominance of the forward scattering processes is the important property which 
yields the correlation function (17) and leads with (15) to the final result (29). 



4 Two channel system 



As yet another application of the general result (15) we now derive the exact 
solution for the case of two open channels N = 2. Here, the matrix u nm is a 



2x2 unitary matrix. Thus only the expression ( 



No 



|wio| occurs in Eq. 



(15) for the localization length. One can rewrite (15) in terms of this variable 
and obtains 



i: 



8D 



- |/?io| 2 )<C 2 ) + 2/3-{C> + + |/?io| 2 ) 



(30) 



with f3 ± = (|/9oo| 2 =t|/3n| 2 )/4. This reduces the ensemble average over combina- 
tions of matrix elements u nm to the moments {() and (( 2 )- In order to find the 
stationary distribution function P(() we use the parameterization ( = cos 29 
and 



u 



( 



\ 



i<t>0 Q 







cos 9 — sin 9 
sin 9 cos 9 




(31) 



for the unitary matrix u nm . The stationary probability distribution can then 
be obtained from the associated Fokker-Planck equation. Analyzing this equa- 
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tion we find that the exponentials e*^ ' 1 and e* 7 always occur together with 
the matrix elements a and /3, respectively. If one neglects the rapidly oscillat- 
ing contributions in this Fokker-Planck equation then only those terms which 
contain the magnitudes |a nm | 2 and |/3 nm | 2 remain. For that reason, the sta- 
tionary probability distribution is independent of the phases 0i 5 2 and 7. The 
solution -P(C) of the stationary Fokker-Planck equation follows from the first 
order differential equation 



I' 1 



with A : 



10 



-AC 2 ; 



B 



P(C) 



|«io| 2 + ^ 



and B 



2/3" 



|«io| 2 + /3+' 



(32) 
(33) 



It is of the type P(() ~ exp[-P J c dx(l + Ax 2 )- 1 }/ {I + AC, 2 ). Thus we arrive 
at the following expression for the two channel localization length: 



SEl-\\a lQ \ 2 + (3' 



1 + 



|/3io| 2 + ^ 
|aio| 2 + /9 J 



«10 


2 + 


P10 


2 




a w 


2 + (3+ 



P+ 



(34) 



with P + = [P(l) + P(— 1)]. The length scale I is the mean free path of a single 
particle having the same energy as the pair, i.e. E = 2Ep. The probabilities 
P(±l) for ( — ±1 can easily be obtained from the solution of Eq. (32). It 
yields 



P + : 
P + 



B 
B 



coth 
coth 



B 
B 

7^ 



arctan VA 



for A > 



artanhV— A 



for A < 0. 



Expression (34) provides a general solution for the localization length in the 
two channel regime without any further restrictions concerning the matrix 
elements a nm and (5 nm . It is therefore the main result regarding the two channel 
case. 

Some special cases of (34) may be of interest. For isotropic backward scat- 
tering \f3 m \ 2 = \p n \ 2 = 2\p 10 \ 2 follows A = and B = in (32) and the 
distribution function is independent of the forward scattering: P(C) = 1/2. 
This case corresponds to the invariant distribution of an ensemble of unitary 
2x2 matrices. If the two channels are equivalent, so that | /?oo 1 2 — \Pn\ 2 an d 
thus B = 0, the distribution function P(() ~ (1 + AC 2 )- 1 follows from (32). 
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-4 -2 ,, 2 4 



Fig. 2. Two channel case N = 2: The two interaction potentials under consideration 
are u osc (y) (thick solid line) and Ub ox (y) (thick dashed line). The width of the box 
and the absolut position of the oscillator potential have been choosen so that the 
energy levels eo,i (thin lines) of the first two symmetric bound states coincide. 

The special case where \/3 00 \ 2 = \(3n\ 2 = \Pw\ 2 was already discussed in [7] and 
is correctly reproduced by our general solution. 

We complete our discussion of the N = 2 case with the presentation of some 
numerical results for two specific examples of interaction potentials. Using the 
result (34) one can analyze the influence of the interaction potential u(y) on 
the localization length l c . The requirement of exactly two open channels re- 
stricts the energy range to e\ < E < €2 where eo, ei, . . . , e n are the energies 
corresponding to the bound states with given parity. In order to demonstrate 
how the localization length depends on the structure of the pair, we present 
a comparison between the harmonic oscillator interaction u osc (y) and the in- 
teraction where the two particles are bound together by a box-like potential 
u box{y)- While u osc (y) is a smooth function of the relative coordinate y the sec- 
ond interaction potential Ub ox (y) is constant for \y\ being smaller than some 
fixed value y s and jumps to infinity outside this region, see Fig. 2. 

We first discuss the case of bound states with even parity. The size of the box 
and a total energy shift for the harmonic oscillator potential were adjusted so 
that the energy levels eo,i entering (34) are at the same position. This approach 
ensures that all direct energy dependencies of the localization length have the 
same structure for either interaction type. The energy window for the N = 2 
case is then given by 4.5 < E£ 2 < 8.5 (4.5 < E£ 2 < 12.5) for the harmonic 
oscillator (box-like potential), Fig. 2. 

The next step is to evaluate the matrix elements a and (5 as given by (12) 
and (13). Inserting these matrix elements in (34) gives the result for the local- 
ization length l c (E) which is presented in Fig. 3. The interaction between the 
constituents of the composed particle clearly causes a derealization compared 
to the single particle with no internal degrees of freedom. The reason is that 
the finite size of the pair effectively smoothes out the disorder. However, this 
derealization effect sensitively depends on the total energy E and the type 
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Fig. 3. Two channel case N = 2: The pair localization length l c compared to the 
single particle localization length I as a function of the dimensionless energy E£ 2 
(even parity of the bound state wave functions 4> n ). 



0.4 



0.2 



osc — 




box - 











Fig. 4. Two channel case N = 2: The ratio between forward scattering length If and 
backward scattering length If, as a function of the dimensionless energy E^ 2 (even 
parity of the wave functions (f) n ). 

of the two particle interaction. The resonance in the case of the box-like po- 
tential is clearly more pronounced. If the energy is close to this resonance, i.e. 
E£ 2 ps 6.6, the derealization effect differs by a factor ~ 5 due to the different 
interaction potentials. 

Furthermore, we have examined the ratio between forward and backward scat- 
tering If /h as a function of the energy. As Fig. 4 indicates, the forward scat- 
tering relative to the backward scattering is already for just two channels 
enhanced. Although the ratio If/h depends on the energy this effect holds 
true for the entire energy window accessible for the two channel case. The 
suppression of the backward scattering is because the finite size of the pair re- 
duces scattering processes with large momentum transfer as can be concluded 
from Eqs. (12, 13). This effect is even more pronounced if many channels are 
open, a fact that is of importance for the solution of the multi-channel case, 
as discussed in Sec. 3. 

Using the sum rule (9) we calculate the length scale l\ of the transmission 
eigenvalue of the second internal channel. It is given by l\ = (21^ — i" 1 ) -1 , 
see Fig. 5. We find that the ratio l\jl c is smaller than 0.6 over the entire energy 
range for either interaction potential. This justifies the assumption that the 
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Fig. 5. Two channel case N = 2: The ration l\/l c of the length scales of the two 
channels as a function of the dimensionless energy E^ 2 (even parity of the wave 
functions (f> n ). 




Fig. 6. Two channel case N = 2: The pair localization length l c compared to the 
single particle localization length I as a function of the dimensionless energy E£ 2 
(odd parity of the wave functions </>«,)• 



transmission is mainly due to one internal channel if the system is large, i.e. 
L>1, because of the exponential dependence of the transmission amplitudes 
on the system size. 

Finally we compute the localization length for the situation where the two 
bound states are of odd parity. Again, the parameters of the interaction po- 
tentials are fixed such that the energies of the two involved energy levels coin- 
cide for the different interaction potentials. The resulting localization length 
is plotted in Fig. 6. There is a significant difference between the localization 
lengths in the case of wave functions with even parity compared to odd parity. 
For the harmonic oscillator potential there is a sharp resonance near the lower 
limit Ei of the energy window in the second case (Fig. 6) while the first case 
shows a smooth behavior with a weak resonance in the middle of the energy 
range (Fig. 3). For the Ub ox (y) interaction the resonance position and height 
also changes significantly. Since the energy levels are shifted to the same po- 
sition for either case the differences result only from the different structure 
of the bound state wave functions. The largest derealization effect occurs for 
u box(y) f° r the states of odd parity and is about l c /l ~ 25 at the resonance. 
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5 Conclusion 



In this work we have studied the localization of a pair of two bound particles in 
a one-dimensional system with weak disorder. The interaction potential within 
this composite particle is fully taken into account within the semiclassical 
approximation. Our results are based on a method introduced by Dorokhov 
for two particles bound by a harmonic oscillator potential. We have generalized 
the method from this special interaction potential to an arbitrary attractive 
interaction. This allows us to analyze the dependence of the pair localization 
on the interaction. We find an enhancement of the pair localization length l c 
in comparison to the single particle mean free path /. This enhancement is 
independent of the form of the pair interaction potential in the limit of many 
bound states N >> 1. It is given by l c /l ~ N as can be seen from the central 
result (29) in conjunction with (24). Furthermore we derived an exact solution 
for the two channel case N = 2. For a bound pair with N — 1, 2 we observe a 
sensitive dependence of the localization length (16), (34) on the shape of the 
interaction potential, the parity of the involved bound states and the kinetic 
energy of the pair. We expect that such a behavior is typical for bound pair 
with a small number of bound states. 

The approach described in this work is valid for two identical particles. Nev- 
ertheless, it is possible to extend the method without difficulties to a pair of 
different particles (e.g. electron-hole pair). In this case, the expression for the 
matrix elements (13) has to be modified taking into account that the random 
potential and the mass may be different for the two particles. However, the 
described method cannot be extended to a repulsing interaction because the 
assumption that the pair size is smaller than the mean free path breaks down. 
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